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Abstract 

A three-bay, space, cantilever truss is 
probabilistically evaluated to describe progres- 
sive buckling and truss collapse in view of the 
numerous uncertainties associated with the 
structural, material and load variables (primi- 
tive variables) that describe the truss. Ini- 
tially, the truss is deterministically analyzed 
for member forces and member(s) in which the 
axial force exceeds the Euler buckling load are 
identified. These member(s) are then discre- 
tized with several intermediate nodes and a 
probabilistic buckling analysis is performed on 
the truss to obtain its probabilistic buckling 
loads and the respective mode shapes. Fur- 
thermore, sensitivities associated with the 
uncertainties in the primitive variables are 
investigated, margin of safety values for the 
truss are determined and truss end node dis- 
placements are noted. These steps are 
repeated by sequentially removing buckled 
member(s) until onset of truss collapse is 
reached. Results show that this procedure 
yields an optimum truss configuration for a 
given loading and for a specified reliability. 

Introduction 

It is customary to evaluate the structural 
integrity of trusses by using deterministic 
analysis techniques and appropriate load/ 
safety factors. Traditionally, these factors are 
an outcome of many years of analytical, as 
well as experimental, experience in the areas 
of structural mechanics/design. Load factors 
are used to take into account for uncertainties 
in many different operating conditions includ- 
ing the maximum loads and safety factors are 
also used to account for unknown effects in 
analysis assumptions, fabrication tolerances, 
and material properties. 


As an alternative to the deterministic 
approach, is the Probabilistic Analysis 
Method (PS AM). 1 This method formally 
accounts for various uncertainties in primitive 
variables (fundamental parameters describing 
the structural problem) and uses different 
distributions such as the Weibull, normal, log- 
normal, etc. to define these uncertainties. 
Furthermore, PSAM assesses the effects of 
these uncertainties on the scatter of structural 
responses (displacements, frequencies, eigen- 
values). Thus, PSAM provide a more realistic 
and systematic way to evaluate structural 
performance and durability. A part of PSAM 
is a computer code NESSUS (Numerical Eval- 
uation of Stochastic Structures Under Stress) 
which provides a choice of solution for static, 
dynamic, buckling, and nonlinear analysis. 2,3 

In the recent past, NESSUS has been 
used for the analysis of Space Shuttle Main 
Engine (SSME) components. Representative 
examples include a probabilistic assessment of 
a mistuned bladed disk assembly 4 and an 
evaluation of the reliability and risk of a tur- 
bine blade under complex service environ- 
ments. 5 Furthermore, NESSUS has also been 
used to computationally simulate and probabi- 
listically evaluate a cantilever truss typical for 
space type structures 6 and quantify the uncer- 
tainties in the structural responses (displace- 
ments, member axial forces, and vibration 
frequencies). The objective of this paper is to 
develop a methodology and to perform proba- 
bilistic progressive buckling assessment of 
space type trusses using the NESSUS com- 
puter code. 

Fundamental Approach and 
Considerations 

One of the major problems encountered 
in the analysis of space type trusses is to come 


up with a stable and optimum configuration 
for given loading conditions and to be able to 
probabilistically analyze them to take into 
account the probable uncertainties in the 
primitive variables typical for space environ- 
ment conditions. Presently, it is a practice 
to design these trusses with cross bracings 
thereby increasing the overall weight of the 
truss, the cost of fabrication and the effort to 
deploy in space. Furthermore, the presently 
available methods/programs do not easily 
allow us to identify any local instability in 
any of the internal members of the truss dur- 
ing probabilistic buckling (eigenvalue) analysis 
and to calculate over all margins of safety of 
the truss. Therefore using the NESSUS code, 
a methodology for the probabilistic progres- 
sive buckling is developed as described 
hereafter. 

Finite Element Model 

A three-dimensional, three-bay cantilever 
truss is computationally simulated using a 


linear isoparametric beam element based on 
the Timoshenko beam equations. The element 
is idealized as a two-noded line segment in 
three-dimensional space. The cantilever truss 
is assumed to be made from hollow circular 
pipe members. The members are made up of 
wrought aluminum alloy (616-w) with modu- 
lus of elasticity (E) equal to 10 Mpsi. The 
outer and inner radii (r and of the tube, 
are 0.5 and 0.4375 in., respectively. All six 
degrees of freedom are restrained at the fixed 
end (left side) nodes. Each bay of the truss is 
5 ft wide, 8 ft long, and 6 ft high (Fig. 1). 

The overall length of the truss is 24 ft. Six 
vertical and two longitudinal loads are 
applied. In addition, twisting moments are 
applied at the truss-end nodes. The directions 
of the forces and moments are shown in Fig. 1 
and mean values are given in Table I. The 
applied loads and moments are selected to 
represent anticipated loading conditions for a 
typical space truss. 


TABLE I. - PRIMITIVE VARIABLES AND UNCERTAINTIES FOR 
PROBABILISTIC STRUCTURAL ANALYSIS OF A SPACE TRUSS 


[Random input data.] 


Primitive variables 

Distribution 

Mean 

Scatter, 



type 

value 

dipercentage 

Geometry 

Width 

Normal 

60 in. 

0.5 


Length 

Normal 

96 in. 

0.1 




192 in. 

.1 




288 in. 

.1 


Height 

Normal 

72 in. 

0.2 

Loads 

Vertical 

Lognormal 

20 lb 

6.3 


Longitudinal 

Lognormal 

20 lb 

2.5 


Twisting moment 

Lognormal 

50 lb in. 

6.3 

Material 

Modulus 

Normal 

10 Mpsi 

7.5 

property 





Tube 

Outer radius 

Normal 

0.5 in. 

7.5 

radii 

Inner radius 

Normal 

0.44 in. 

7.5 
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Buckling of Columns 

In slender columns, a relatively small 
increase in the axial compressive forces will 
result only in axial shortening of the member. 
However, the member suddenly bows out side- 
ways if the load level reaches to a certain 
critical level. Large deformations caused by 
increased induced bending moment levels may 
lead to the collapse of the member. On the 
other hand, tension members as well as short 
stocky columns fail when the stress in the 
member reaches a certain limiting strength of 
the material. According to Chajes 7 "...Buck- 
ling, however, does not occur as a result of the 
applied stress reaching a certain predictable 
strength of the material. Instead, the stress 
at which buckling occurs depends on a variety 
of factors, including the dimensions of the 
member, the way in which the member is sup- 
ported, and the properties of the material out 
of which the member is made...". Chajes also 
describes the concept of neutral equilibrium 
which is being used to determine the critical 
load of a member such that at this load level 
the member can be in equilibrium both in the 
straight and in a slightly bent configuration. 
Furthermore, the Euler load (buckling load or 
critical load) is the smallest load at which a 
state of neutral equilibrium is possible or the 
member ceases to be in stable configuration. 
This above definition of buckling load is used 
to identify the probable truss members that 
contribute to the progressive buckling behav- 
ior of the cantilever truss. 

Probabilistic Model 

The following primitive variables are 
considered in the probabilistic analysis: 

(1) nodal coordinates (X,Y,Z) 

(2) modulus of elasticity (E) 

(3) outer radius of the tube (r ) 

(4) inner radius of the tube (r^ 

(5) vertical loads (V) 


(6) longitudinal loads (H) 

(7) twisting moments (M) 

It is possible that the above primitive 
variables will vary continuously and simulta- 
neously due to extreme changes in the envi- 
ronment when such trusses are used in upper 
Earth orbit for space station type structures. 
The normal distribution is used to represent 
the uncertainties in E, r Q , r i? and X,Y,Z 
coordinates. The applied loads and moments 
are selected to represent an anticipated load- 
ing for a typical space truss. The scatter in 
these are represented by log-normal distribu- 
tions. Initially, the NESSUS/FEM (Finite 
Element Methods) module is used to deter- 
ministically analyze the truss for mean values 
of each of these primitive variables. In the 
subsequent probabilistic analysis, each primi- 
tive variable is perturbed independently and 
by a different amount. Usually, the perturbed 
value of the primitive variable is obtained by 
certain factor of the standard deviation on 
either side of the mean value. It is important 
to note that, in the NESSUS code a linear 
buckling analysis is carried out by making use 
of the subspace iteration technique to evaluate 
the probabilistic buckling load. The matrix 
equation for the buckling (eigenvalue) analysis 
for a linear elastic structure is as follows: 

{[K] - A[KJ}{^> _ 0 « 

In the above equation, [K] is the standard 
stiffness matrix, [K g ] is the geometric stiffness 
matrix, A is the eigenvalue, and <f> are the 
eigenvectors. Finally, the NESSUS/FPI (Fast 
Probability Integration) module extracts 
eigenvalues to calculate a probability distri- 
bution of the eigenvalues and to evaluate 
respective sensitivities associated with the 
corresponding uncertainties in the primitive 
variables. The mean, distribution type and 
percentage variation for each of the primitive 
variables are given in Table I. 
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Probabilistic Progressive Buckling 
Computational Simulation 

Initially, the truss is deterministically 
analyzed for member forces and identify the 
member(s) in which the axial forces exceed 
the Euler load. These member(s) were then 
discretized with several intermediate nodes 
and a probabilistic buckling (eigenvalue) anal- 
ysis is performed to obtain probabilistic buck- 
ling loads and respective buckled shapes. 
Furthermore, the sensitivity factors represent- 
ing the impact of uncertainties in the prim- 
itive variables on the scatter of response 
variable (eigenvalue) are evaluated. Finally, 
any member(s) that have buckled are identi- 
fied and the probabilistic buckled loads/ 
moments at each probability level are 
obtained by multiplying the respective eigen- 
values with the applied loads and moments. 

In the subsequent analyses the buckled 
member(s) are removed from the original truss 
configuration and the above described analysis 
steps are repeated until onset of collapse state 
is reached. It is important to note that the 
mean values of the loads and moments are 
kept constant and are perturbed around their 
means during the probabilistic buckling analy- 
sis. The truss end node displacements versus 
the number of members removed are plotted 
to identify the onset of the truss collapse 
state. Finally, the minimum number of mem- 
bers needed to support the applied loads and 
moments are determined. 

Discussion of Results 

Probabilistic Progressive Buckling - First 
Buckled Member 

Figures 2(a) to (f) show the probabilistic 
progressive buckled mode shapes of the three- 
bay space truss as individual buckled members 
are sequentially removed from the original 
configuration until it reaches the onset of 
collapse. The probabilistic buckling analysis 
indicated that the first bay front diagonal 
buckled first, (Fig. 2(b)) and the correspond- 
ing probabilistic buckled loads and moments 


at 0.5 probability are shown for example in 
Fig. 3. Probabilistic buckled loads and mo- 
ments at different probability levels can also 
be obtained. Furthermore, a method of calcu- 
lating the margin of safety (MOS) for speci- 
fied probability by using known distributions 
for applied loads and moments and corre- 
sponding cumulative distribution function 
curves obtained from PS AM are shown in 
Fig. 4. The sensitivity factors from Fig. 5 
suggest that the scatter in the bay length 
parameter (Y-coordinate) had the highest 
impact on the probabilistic distribution of the 
buckling load followed by the bay height 
(Z-coordinate), bay width (X-coordinate), 
vertical and longitudinal loads and finally 
twisting moments. Any slight variation in 
spacial (geometry) variables has a direct effect 
on the overall length of the members and 
thereby alters many terms in the stiffness 
matrix containing the length parameter. 
Finally, this has a definite affect on the proba- 
bilistic buckling loads which has been clearly 
observed in the above discussed results. How- 
ever, it is important to note that even com- 
paratively large variations in both member 
modulus (E) and area (r Q and rj) (see 
Table I) had very negligible impact. Similar 
conclusions can also be drawn for the probabi- 
listic member force in the first buckled 
member (see Fig. 6). The variation in the 
resistance (mean area x mean yield strength) 
of the member was assumed to have a Weibull 
distribution and is shown in Fig. 7. MOS 
calculations for strength exceedence using 
distribution curves for probabilistic member 
force and resistance as well as probabilistic 
buckling load and resistance indicate that the 
buckled member did satisfy the strength cri- 
teria condition. Therefore, it can be con- 
cluded from Figs. 4 and 7 that the member 
buckled when its axial force exceeded the 
Euler buckling load and when the stress due 
to this load did not exceed the failure criteria. 

Probabilistic Progressive Buckling - 
Second/Third/Fourth Buckled Members 

As described in the previous section, the 
deterministic analysis followed by the 


4 



probabilistic analysis was performed with 
sequential removal of the first, second, third, 
and fourth buckled member from the truss 
and the probabilistic buckled loads and 
moments, sensitivity factors, and MOS values 
for stress were obtained. When the second 
member wets buckled (see Fig. 2(c)) the com- 
parable results as shown in Figs. 4 to 7 are 
described in Figs. 8 to 11. For these truss 
configurations the MOS value decreased 
from 3.53 to 2.53. The similar details of the 
truss with the third buckled (see Fig. 2(d)) 
are given in Figs. 12 to 15. According to 
Fig. 13, the MOS value further decreased to 
1.62. Similarly, Figs. 16 to 19 give compara- 
ble results of the truss when the fourth mem- 
ber was buckled (see Fig. 2(e)). It is 
important to note from Figs. 18 and 19 that, 
the scatter in the bay height had much higher 
impact than scatter in both bay width and 
length on both the probabilistic buckling 
loads/mornents and buckled member force. 
Finally, the details of the onset of collapse 
state of the truss (Fig. 2(f)) are shown in 
Figs. 20 to 22. When all the four buckled 
members were removed, the MOS value was 
equal to -3.75 which indicates that the onset 
of collapse was reached (Fig. 21). Further- 
more, the probabilistic buckling loads/ 
moments at 0,001 probability were equal to 
maximum applied loads/moments with 
assumed distributions (see Fig. 20). In addi- 
tion, at the collapse state the uncertainties in 
both the bay length and bay height had suffi- 
ciently high impact on the probabilistic buck- 
ling loads/moment distributions (see Fig. 22). 
In the above discussed various truss configura- 
tions, the uncertainties in the vertical loads 
had consistently the same impact on buckling 
loads/moment, where as member modulus and 
area had negligible impact. 

Probabilistic Truss End Node Displacements 

The truss end node displacements (later- 
al, longitudinal, and lateral) were also cal- 
culated during each of above mentioned 
deterministic analyses for each truss configu- 
ration and are shown in Fig. 23. It is clear 
that there is not considerable change in either 


lateral or longitudinal displacement as each 
buckled member was sequentially removed. 
However, the truss end node vertical displace- 
ment gradually increased up to the truss con- 
figuration with three buckled members 
removed and suddenly increased very rapidly 
when the fourth buckled member was removed 
giving an indication of unbounded displace- 
ment growth which suggests that the truss 
had reached the onset of its collapse state. 

This is due to the fact that the total vertical 
loads are six times higher than total longitudi- 
nal loads and the perturbations in the vertical 
loads are higher than that of twisting 
moments. Figures 24 and 25, respectively, 
show the relationships between the applied 
vertical loads and probabilistic buckling loads 
as well as probabilistic buckling loads and 
MOS values. The optimum truss configu- 
ration was reached with the forth buckled 
member removed whereby the probabilistic 
buckling load was equal to the applied vertical 
load at 0.001 probability level (see Fig. 24). 
Similar conclusions can also be made for longi- 
tudinal loads and twisting moments. In addi- 
tion, there is a gradual decrease in the MOS 
values as buckled members were sequentially 
removed and reached a zero value when the 
optimum truss configuration was reached (see 
Fig. 25). Similar conclusions can also be 
made for longitudinal loads and twisting 
moments. 

Probabilistic Buckling Including Initial 
Eccentricity 

In the above discussed probabilistic pro- 
gressive buckling methodology, all the mem- 
bers were assumed to be initially perfectly 
straight and the buckled members were 
sequentially removed with the assumption 
that once the member buckled it would yield 
and could not resist any additional loading 
and thereby would not contribute to the over- 
all stiffness of the truss. In order to verify 
this assumption, the maximum eccentricity at 
which the yielding in the member (first bay 
front diagonal) will take place due to the 
combined effects of axial and in-plane bending 
moments was calculated. Furthermore, this 


5 


member was modeled to depict the buckled 
configuration of the member at which yielding 
will take place, using a parabolic distribution 
for the above calculated eccentricity (see 
Fig. 26). The deterministic and subsequent 
probabilistic buckling analyses indicate, 
respectively, that the probabilistic buckling 
loads and moments did not change signifi- 
cantly from the original analysis (see Fig. 3) 
and the first bay rear diagonal has buckled 
(see Fig. 27). However, as seen from Figs. 5 
and 28 for probabilistic buckling loads and 
from Figs. 6 and 29 for probabilistic member 
forces, the sensitivity factors show some 
changes especially the variations in bay width 
has the most dominant impact on both proba- 
bilistic buckling loads and moments (see 
Figs. 28 and 29). This is due the fact that 
member buckles in the plane perpendicular to 
the direction of the loading. Nevertheless, it 
is important to note that the scatter in the 
spacial location accentuates the sensitivities of 
the bay length/width/height on the probabil- 
istic load and diminishes that of vertical load. 
Once again the variations in the member 
modulus and area have very negligible impact. 
These results justify the sequential removal of 
the buckled members during progressive 
buckling. 

Conclusions 

The computational simulation of probabi- 
listic progressive buckling of a space type 
truss is demonstrated using the NESSUS com- 
puter code and step-by-step procedure is out- 
lined. The methodology for onset of truss 
collapse is established. The progressive buck- 
led mode shapes are obtained and sensitivities 
associated with the uncertainties in the primi- 
tive variables are evaluated. The safety mar- 
gins are determined and the optimum truss is 
obtained for the given loading condition. The 
results indicate that: (1) probabilistic buck- 
ling loads and margin of safety values decrease 
as buckled member (s) are sequentially re- 
moved; (2) the scatter in truss geometry (bay 
length/wid th/height) and vertical loads have 
considerable impact on the probability of the 


the buckled load; (3) the member modulus 
and area parameters have negligible impact; 
and (4) initial eccentrics have negligible influ- 
ence on the probabilistic buckling load but 
may influence the sensitivities. Collectively 
the results demonstrate that the probability of 
collapse of space-type trusses can be reliably 
assessed by the procedure described herein. 
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(C). (f). 

Figure 2. — Probabilistic progressive buckling as buckled members sequentially removed. 


First buckled member - First bay front diagonal 



Probabilist ic buckled loading condition (Q.5 Prob.) 

Vertical load 782 lb 
Longitudinal load 260 lb 
Twisting moment 651 lb-in. 

Figure 3.-— Probabilistic progressive buckling. 
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Cumulative probability 



Vertical 

load 

Longitudinal 

load 

Twisting 

„ moment 

Applied Response 



Vertical/Longitudinal load, lb. 
Twisting moment, Ib.-in. 



Vertical load 782 lb 
Longitudinal load 2601b 
Twisting moment 651 lb-in. 


* 

5 


£ 


1.00 i— 

.75 — 



Bay length 
Bay width 
Bay height 
Vertical load 
Other load 

Member modulus and 
area negligible 



Probability level for buckling 


Figure 4. — Probabilistic buckling load. 


Figure 5.— Probabilistic buckling load sensitivity. 



Vertical load 782 lb 
Longitudinal load 260 lb 
Twisting moment 651 lb-in. 


BBS Bay length 
i i Bay width 
r~i Bay height 
Vertical load 



0.0238 0.9999 


Probability level for buckling 





Member force/ 
Resistance, lb 


Buckling load/ 
Resistance.lb 


MOS « R at 0.001 -1 
PMF at 
0.999 

MOS ■ 3000-1 - 24.4 
123 


MOS - R at 0.001 -i 
PBLat 
0.999 


MOS « 3000-1 « 3.75 
801 


Figure 6. — Probabilistic member force sensitivity. 


Figure 7. — Probability of strength exceedence (for first bay 
front diagonal). 
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First buckled member removed - First bay front diagonal 

Second bucled member - First bay rear diagonal 



Vertical load 566 lb 
Longitudinal load 1 96 lb 
Twisting moment 486 lb-in. 



Figure 8. — Probabilistic progressive buckling. 


Figure 9. — Probabilistic buckling load. 



Vertical load 566 lb 
Longitudinal load 1 96 lb 
Twisting moment 488 lb-in. 



Vertical load 586 lb 
Longitudinal load 1 96 lb 
Twisting moment 468 Ib-ln. 



0.0241 0.9999 


Probability level for buckling 
Figure 10. — Probabilistic buckling load sensitivity. 



0.0241 0.9999 

Probability level for buckling 


Figure 1 1 . — Probabilistic member force sensitivity. 
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Sensitivity 


First buckled member removed - First bay front diagonal 
Second buckled member removed - First bay rear diagonal 
Third buckled member - Second bay rear diagonal 



Vertical load 440 lb 
Longitudinal load 147 lb 
Twisting moment 367 lb-in. 


Figure 1 2.— Probabilistic progressive buckling. 



Vertical load 440 lb 
Longitudinal load 1471b 
Twisting moment 367 lb-in. 


BBSS Bay length 



0.0239 0.9999 

Probability level for buckling 


Figure 14. — Probabilistic buckling load sensitivity. 
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Sensitivity Cumulative probability 



Applied Response 



0 250 500 750 1 000 

Verticai/Longitudinai load, lb. 

Twisting moment, lb-in. 


Figure 13.— Probabilistic buckling load. 



Vertical load 440 lb 
Longitudinal load 147 lb 
Twisting moment 367 Ib-ln. 



0.0239 0.9999 

Probability level for buckling 


Figure 1 5.— Probabilistic member force sensitivity. 




First buckled member removed - First bay front diagonal 
Second buckled member removed -First bay rear diagonal 
Third buckled member removed - Second bay rear diagonal 
Fourth buckled member - Second bay front diagonal 



Probabil istic buckled loading condition 1Q.5 Prob.I 


Vertical load 233 lb 
Longitudinal load 78 lb 
Twisting moment 1 94 lb-in. 



Applied Response 



Vertical/Longitudinal load, lb. 
Twisting moment, Ib.-in. 


Figure 16. — Probabilistic progressive buckling. 


Figure 1 7.— Probabilistic buckling load. 



Vertical load 233 lb 
Longitudinal load 78 lb 
Twisting moment 194lb-ln. 



Vertical load 233 lb 
Longitudinal load 78 ib 
Twisting moment 194 Ib-ln. 


BBS Bay length 
Bay width 



1.00 


.75 



.25 

0 



Bay length 
Bay width 
Bay height 
Vertical load 
Other load 



0.0247 


0.9999 


Probability level for buckling 


Figure 18. — Probabilistic buckling load sensitivity. 


Figure 19.— Probabilistic member force sensitivity. 



First buckled member removed - First bay front diagonal 
Second buckled member removed - First bay rear diagonal 
Third buckled member removed - Second bay rear diagonal 
Fourth buckled member removed - Second bay front diagonal 

Onset of collapse state 



Vertical load 1401b 
Longitudinal load 47 lb 
T wisting moment 117 Ib-ln. 


Figure 20. — Probabilistic progressive buckling. 


Vertical load 1401b 
Longitudinal load 47 ib 
Twisting moment 117 lb-in. 



0.0227 0.9999 

Probability level for buckling 


Figure 22. — Probabilistic buckling load sensitivity. 
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Cumulative probability 



Applied Response 



Vertical/Longitudinal load, lb. 
Twisting moment, Ib.-ln. 


Figure 21 .—Probabilistic buckling load. 



Number of members removed 


Figure 23.— Progressive buckling leading to structural collapse 
as indicated by unbounded displacement. 




Figure 24.— Buckling load vs applied load/buckiing load at 0.001 
probability. 



Figure 26.— First bay front diagonal with initial eccentricity. 



Figure 25. — Buckling load vs margin of safety at 0.001 
probability. 



Probabilistic bucked loading condition (0.5 prob.) 

Vertical load 762 lb 

Longitudinal load 254 lb 

Twisting moment 635 lb-in. 

Figure 27. — Probabilistic buckling-flrst bay rear diagonal buckled. 



Vertical load 762 lb 
Longitudinal load 254 lb 
Twisting moment 635 Ib-ln. 



Vertical load 762 lb 
Longitudinal load 254 lb 
Twisting moment 635 Ib-ln. 


1.00 , — 


> 

.ti 

(0 


-75 \— 


.50 


-25 \— 


Bay length 
Bay width 
Bay height 
Vertical load 
Other load 

Member modulus and 



0.0238 0.9999 

Probability level for buckling 

Figure 28. — Probabilistic buckling load sensitivity. 


Bay length 
Bay width 



Figure 29. — Probabilistic member force sensitivity. 
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